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fields  are  increasing  monotonically  with  time.  The  general  character  of  the 
near* tip  fields  Is  analyzed  both  with  respect  to  its  variation  with  time  and 
with  polar  angle.  The  non-linear  near-tip  fields  are  related  to  the  linearly 
elastic  far-fleld  by  means  of  a path- Independent  Integral.  In  the  second  part 
of  the  paper  we  consider  rapidly  propagating  cracks.  We  discuss  the  near- 
tip fields  for  various  models  of  material  behavior.  In  particularly  we 
briefly  review  some  earlier  work  by  Achenbach  and  Kanninen  for  a rapidly 
propagating  Mode-Ill  crack,  in  a material  which  displays  strain  hardening. 

In  the  last  part  of  the  paper  we  consider  the  fields  near  a rapidly  propa- 
gating crack-tip  in  an  elas tic-perfect ly  plastic  material  for  the  case  that 
inertial  terms  are  of  importance.  The  system  of  governing  equations  in  the 
plastic  region  is  presented  and  shown  to  be  hyperbolic  in  nature.  As  a first 
approximation  the  steady-state  case  with  respect  to  the  moving  crack-tip  is 
considered  and  an  asymptotic  analysis  of  the  near-tip  field  i3  carried  out. 
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NEAR- TIP  PLASTIC  DEFORMATIONS  IN  DYNAMIC  FRACTURE  PROBLEMS 


J.  D.  Achenbech,  P.  Burger i and  V.  Dunayevsky 
Tha  Technological  Institute 
Northwestern  University 
Evanston,  111.  60201 


ABSTRACT 

Under  rapid  loading  conditions  and/or  for  a rapidly  propagacing  crack, 
tha  mass  density  of  a material  affects  the  fields  of  stress  and  deformation. 
For  such  dynamic  fracture  problems  plastic  deformations  in  the  icmedlate 
vicinity  of  a crack  tip  are  Investigated  In  this  paper.  Both  stationary 
and  propagating  crack  tips  are  considered.  For  a stationary  crack  tip, 
detonation  theory  la  employed  for  the  first  phase  of  the  loading  whan  the 
flelda  are  Increealng  monotonlcally  with  time.  The  general  charactar  of  the 
near- tip  fields  Is  analysed  both  with  respect  to  Its  variation  with  time  and 
with  polar  angla.  The  non-llneer  near- tip  fields  are  relaced  to  the  linearly 
elastic  far* field  by  means  of  a path -independent  integral . In  tha  second  part 
of  the  paper  ve  consider  rapidly  propagating  cracks.  He  discuss  the  near- 
tip fields  for  various  models  of  material  behavior.  In  particularly  ve 
briefly  review  some  earlier  work  by  Achanbach  and  Kanninan  for  a rapidly 
propagating  Mode-Ill  crack.  In  a material  which  displays  strain  hardening. 

In  tha  last  part  of  the  paper  ve  consider  the  fields  near  a rapidly  propa- 
gating crack-tip  in  an  elastlc-perfectly  plastic  material  for  the  case  that 
Inertial  terms  are  of  Importance.  The  system  of  governing  equations  in  the 
plastic  region  is  presented  and  shown  to  be  hyperbolic  in  nature.  As  a first 
approximation  the  steady  state  case  with  respect  to  the  moving  crack-tip  is 
considered  and  an  asymptotic  analysis  of  the  near- tip  field  is  carried  out. 

INTRODUCTION 

There  are  two  broad  classes  of  fracture  mechanics  problems  that  may  have 
to  be  treated  as  dynamic  problems.  These  are  concerned  with:  (1)  cracked 
bodies  subjected  to  rapidly  varying  loads.  (2)  bodies  containing  rapidly 
propagating  cracks.  In  both  cases  the  crack  dp  la  in  an  environment  of 
rapidly  varying  fialds  of  stress  and  deformation. 

Impact  and  vibration  problems  fall  Into  the  first  class  of  dynamic 
problems.  In  the  analysis  of  such  problams  it  is  oftan  found  that  the 
dynamic  stresses  near  flaws  are  higher  than  the  stresses  computed  from  the 
corresponding  problem  of  static  equilibrium,  the  dynamic  stress  "overshoot" 
can  be  especially  pronounced  for  cracks.  In  view  of  the  dynamic  amplifica- 
tion, It  Is  conceivable  that  there  are  cases  for  which  fracture  at  a crack 
dp  does  not  occur  under  a gradually  applied  system  of  loeds,  but  where  a 
crack  does  Indeed  propagate  when  the  same  system  of  loads  Is  rapidly  applied. 
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and  gives  rise  Co  waves,  which  strike  Cha  crack  dp. 

Tha  second  class  of  problems  is  squally  important,  since  thera  are 
several  kinds  of  large  engineering  structures  in  which  rapid  crack  growth  i« 
a definite  possibility.  Whan  a crack  propagates  rapidly,  dynamic  effects 
affact  cha  stress  fields  near  cha  crack  tip,  and  henca  tha  conditions  for 
further  unstable  crack  propagation  or  for  crack  arrest. 

In  recant  yaars  there  have  baan  a number  of  comprehensive  review  articles 
in  Che  general  area  of  elastodynamlc  fracture  mechanics,  see  Refs . [l]-L5 ] . 

At  present,  dynamic  fracture  mechanics  solutions  are,  however,  largely  con- 
fined to  conditions  where  linear  elastic  fracture  mechanics  (LEFM)  is  valid. 

The  a las tic -plastic  treatments  required  to  give  accurate  results  have  not  yet 
been  developed  in  a completely  acceptable  manner,  even  under  static  conditions. 
Current  progress  in  this  area,  and  a starting  point  for  the  development  of  a 
dynamic  plastic  propagating  crack  tip  analysis  have  recently  been  discussed 
in  Ref. [6]. 

In  this  paper  further  investigations  are  reported  on  plastic  deformations 
in  the  immediate  vicinity  of  a crack  tip  for  both  stationary  and  moving  crack 
tips. 

The  analysis  for  a stationary  crack  tip  is  based  on  deformation  theory. 

For  tha  corresponding  class  of  static  problems  deformation  theory  was  applied 
by  Hutchinson  L2l  and  Rice  and  Rosengren  l8].  In  Section  3 some  of  the 
results  of  [2l  and  !_8]  are  extended  to  dynamic  near- tip  fields.  Deformation 
theory  is  not  valid  for  unloading,  and  consequently  the  results  generally 
apply  only  in  an  initial  time  interval  when  the  fields  of  stress  and  deforma- 
tion Increase  monotonies lly.  Also,  since  tha  nonlinear  near-tip  fields  are 
related  to  the  linear  far-field  by  means  of  a path- independent  Integral,  the 
plastic  deformation  must  be  confined  to  a small  tone  near  tha  crack  tip. 

In  Sections  4 and  3 plastic  deformations  near  a rapidly  propagating 
crack  tip  are  considered,  for  Che  special  case  of  Mode  III  deformations.  The 
analysis  cakes  tha  inertia  terms  lnco  account,  but  it  is  assumed  that  Che 
fields  are  steady-state,  i.e.,  they  are  constant  with  respect  to  an  observer 
who  moves  with  the  crack-tip.  The  corresponding  quasi-static  problem  has  been 
considered  by  Chltaley  and  McClintock  [9]  for  an  elastic  perfectly -plastic 
material,  and  by  Amaelgo  and  Hutchinson  .10]  for  the  case  of  J2-flow  theory 
with  a bilinear  effective  stress-strain  curve. 

Section  4 is  primarily  concerned  with  a discussion  of  the  influence  of 
various  constitutive  behaviors  on  the  neer-tip  fields.  In  this  section  we 
also  briefly  review  some  earlier  results  of  Achenbach  and  Ranninen  [6j  for  32- 
flow  theory.  In  particularly  we  investigate  the  nature  of  the  governing 
equations  as  the  crack-tip  velocity  Increases.  For  the  case  of  strain  hard- 
ening the  governing  equations  are  elliptic  when  the  crack-tip  velocity  is  not 
too  large.  In  chat  case  the  usual  separation  of  variables  asymptotic  analy- 
sis yielding  singular  stress  and  strain  fields  can  be  carried  out.  and 
singularities  of  the  type  rP(-  1 < p < 0)  are  obtained.  As  the  velocicy 
Increases  (or  Che  strain-hardening  curve  becomes  flatter)  the  nature  of  Che 
equations  becomes  hyperbolic.  It  appears  very  difficult  to  trace  the  transi- 
tion of  the  near-tip  fields  as  Che  velocity  Increases. 

UNEAR  E LAS TCDYNAMIC  STRESS  INTENSITY  FACTORS 

It  is  well  known  that  for  linearly  elastic  materials,  Che  stress  fields 
in  the  immediate  vicinity  of  a stationary  or  moving  crack  cip  can  be  ex- 
pressed in  the  general  form 

k(t,v)  Zif(8,v)  , (1) 

where  r,S  are  polar  coordinates  centered  at  the  crack  tip.  as  shown  in  Fig.l. 
v is  the  instantaneous  velocity  of  the  crack  tip  and  is  the  Cauchy 
stress  tensor.  The  functions  I,,(5,v)  are  universal  functions  in  chac  chey 
are  independent  of  the  overall  geometry  of  the  body  and  particular  loading 
systems.  It  is  only  through  a single  parameter,  the  stress  intensity  fector 
k(c,v),  chac  che  overall  geometry  and  loading  of  the  body  Influence  the  near 
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dp  5 crass  fields.  The  form  given  in  Eqn.(l)  is  valid  for  propagating  cracks 
and  stationary  cracks  (v  ■ 0)  under  both  dynamic  and  quasi-static  loading. 

It  is  also  applicable  to  three  dimensional  problems  of  a plane  crack  with  a 
smoothly  curved  edge,  if  the(r,8)-plana  is  taken  normal  to  the  crack  edge. 

The  radius  of  curvature  of  the  crack  edge  then  entera  of  course  in  the 
magnitudes  of  the  stress  intensity  factors.  The  components  of  £.  (5,v)  show 

a significant  dependence  on  v. 

Ic  la  of  note  that  for  the  Mode- 
I case  the  maximum  of  £ga(8,v) 
y i r moves  out  of  the  plane  of  the 

2 T / crack  uhen  v increases  beyond  a 

I / certain  value. 


1 1 Fig.  1 Coordinate  system  for 

a propagating  crack  tip 

DYNAMIC  FIELDS  FOR  STATIONARY  CRACK  TIPS 
USING  DEFORMATION  THEORY 

In  this  section,  an  lncompresalble  elastic  power-law  hardening  material 
is  considered.  The  analysis  is  carried  out  for  the  case  of  plane  strein  and 
the  usual  small  strain-displacement  relations  are  used. 

The  dynamic  Mode-I  stress  Intensity  factor,  k.(t),  for  a stationary 
crack  tip  in  a linearly  elastic  material  can  be  expressed  as 

kT(t)  - f(t)  Kt  (2) 

where  K_  is  the  corresponding  quasi-static  stress  intensity  factor.  It  will 
now  be  shown  that  the  dynamic  screes  intensity  factor  for  the  non-linear 
material  can  be  expressed  in  terms  of  a function  of  f(c),  which  Is  defined 
by  Eq. (2) , and  the  corresponding  non-linear  static  fields. 

Near  Tip  Fields 

Following  Ref.i8],  an  incompressible  material  is  considered,  with  a 
relation  between  devletorlc  stresses,  s^  and  strains, of  the  form 

stJ  - <2t/V)  €tj  (3) 

where  r » (s.,s../2)^  and  y ■ (2*.  c)^.  The  hardening  is  assumed  to  be 
governed  by 

t ■ n v - ( - ~)  y for  y s v 


for  y * yo  , 


where  t , v are,  respectively,  the  yield  stress  and  strain  in  shear,  and  N 
is  the  ftardSning  exponent  satisfying  0 < B s 1.  Since  the  plastic  strains 
are  incompressible  and  only  the  dominant  singularity  Is  considered,  the 
assumption  of  material  incompressibility  should  not  be  a poor  approximation 
for  moderate  strain  hardening.  The  system  of  governing  equations  is  com- 
pleted with  the  equations  of  motion 


To  obtain  eh#  dominating  singularity  at  tha  crack  tip,  tha  naar  tip 
flald  la  axsended  In  earma  of  powers  of  r and  functions  of  c and  9.  The  dis- 


placements ara  expressed  as 

ut  - k4(t)  ut(r,9)  - kt(t)  0^(3)  rq  as  r - 0 (6) 

where  q > 0.  After  use  of  Eqns.(3)  and  (4),  corresponding  expressions  for 
the  strasaas  and  strains,  ara  obtained  as 

V)  ~kc(e)  ®ij(t-9)  - k,(t)  Etj(8)  r(Q*l)N  , (7) 

*lj  " kc(e)  *ij(r'9>  ’ V*)  6^(6)  f,“l  • (8> 

as  r - 0 , where  k (t)  Is  the  "dynamic  stress  Intensity  factor"  and  k (t), 
k((t)  are  related  By 

V°  ■ wlc.Cc)  JM.  (9) 


Substituting  Eqns.(6)  and  (9)  Into  the  aquations  of  motion  (S)  shows 
that  tha  highest  order  singularities  on  the  left  hand-side  of  Eqn.(S)  are  of 

order  vhti*  those  on  the  right-hand  side  are  of  order  rq.  This 

implies,  as  in  the  linear  elastic  case,  that  for  a stationary  crack  the 
inertia  tanas  do  not  affect  tha  variation  with  respect  to  r and  3 of  the 
near- tip  fields  to  the  highest  order  of  singularity.  The  results  for  q, 
u«(9),  £ij(9)  and  E^a(3)  ara  the  same  for  static  and  dynamic  problems,  i.e. 
the  3.R.R.  singular  fields 

*tJ  • 0(r-8^  and  ^ - 0(r^1+K)), 

dominate  in  the  vicinity  of  tha  crack  tip. 

Specific  results  for  the  functions  u.(8),  £. .(9)  and  E..(8)  are  given  in 
Refs. [7]  and  [8].  1 ij  lj 


A Path- Independent  Integral 

For  the  case  of  linear  materials,  Sllsson  [ll] 
independent  Integral,  of  the  fora 

li  “ Jgdl  '’lj*! J * I ; p “iui)nl  ' *1  Ix^]  ds 


has  presented  a path- 


(10) 


in  terms  of  Laplace  transforms  of  tha  field  variables.  The  contour  C defines 
a path  enclosing  the  crack  dp,  n.  ara  the  components  of  the  outer  normal  to 
C,  a superscript  bar  indicates  the  Laplace  transform  and  p is  the  Laplace 
transform  variable. 

For  material  behavior  described  by  nonlinear  deformation  theory,  and 
in  a region  where  Eqns.(6),(7)  and  (8)  ere  valid,  a functional  W(c  ) can  be 
defined  by  * 

W(V’CJ  (li: 

The  following  integral 

XnX  " Jc  [(ks  kt  W<‘lj>  + I 0 pZ  “l  “i>  nl  V*1  ] is  (l2 

ls  chen  path  Independent  as  r ~ 0.  To  prove  this  we  consider  a closed  con- 
tour C , with  enclosed  area  A*  and  we  let  I represent  1^  along  C . 


\ 
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By  eh*  divergence  cheorem, 

-*  r fr  r ^‘ljj 


I.-D 


lc  k 
KC  c 


Jx, 


+ 1 p p 


2 

5x, 


<ui“i> 


A (^ij  3]  «**  (13) 


-J  ' lJ 


where  * 7 Oj  he*  been  used. 


Now 


_ _ »(*-,)  _ _ 3*t1  _ a 

k k ■ ■ j v ■ Ji  ■ j — 

“e  “c  Jxj  Jtj  *«  3«j  ’ij 


j _£i) 

;*i  «y 


* P4,  - 2«(‘km)/:«t . end  the  symmecry  of  3^ 
of  Sqn.(l4)  Into  Eqn.(13)  end  us*  of  Che  Laj 


hes  been  used . 


(14) 

Substitu- 


where 

fcl  M*»  r *>  J 

elon  of  Eqn.(l4)  Into  Eqn.(^3)  end  us*  of  Che  Leplec*  cransform  of  ch* 
equeclon  of  motion,  shows  I • 0 . Reelislng  chec  ch*  fields  used  apply  only 
as  r - 0 , chis  implies  Chat  I , is  path  independent  in  the  vicinity  of  ch* 
crack  tip.  “ 

The  assumption  of  small  scale  yielding  is  now  introduced,  l.e.  ch* 
region  of  nonlinear  deformation  is  assumed  small  as  compared  to  the  region  in 
which  the  linear  elastic  singularity  dominates,  which  in  turn  is  small  with 
respect  to  any  characteristic  dimension  of  the  body  [_12] . Since  I , is  valid 
in  both  the  linear  and  nonlinear  regions  close  to  the  crack  tip,  it  can  be 
used  to  connect  the  nonlinear  elastodynamic  near- tip  field  to  the  surrounding 
linear  elastodynamic  field. 

For  the  linear  elastic  case,  I,  can  be  expressed  in  terms  of  Che  Laplace 
cransform  of  the  dynamic  scress  intensity  factor  ill!.  For  plan*  strain: 


Pl<P>]2 


(1  - v )/E  , 


(15) 


and 


ltu  ’ Vp)  vp)  3 + p2  [vp)l 


2 L 


(16) 


where 


hi. 


and  T • 
vanishes 


Jc[w(‘u)  ai*?iS^]d-’ 

1 " 2 ° I„  “l  “l  *2 

In  the  limit  as  C shrinks  on  to  the  crack  tip,  L 


(17) 

(18) 


C 

k (p)  T, 


Intensity  Faccors 


By  ch*  reasoning  given  above,  w*  have  I,  * I . which  gives 

l tu 


Up)  *f(p)  J • [^(p)]2  (i  - ^2)  / e • 


(19) 


Anologously  co  Eqn.(2)  w*  Introduce 

{k.(p);  kt(p)  } - { fc(p)  ks  ; ft(p)  kf  } 


(20) 


where  K_  , K#  are  the  intensity  factors  for  the  corresponding  scatlc  problem. 

- v‘)/  E and  Eons 


Using  tfi*  quisl-statle  result  that  K K "J  « K?  (1 
(18).  Eqn.(19)  is  reduced  to 


f.(p)  ft(P>  - [f(p)]2 


Eons  (2  "i  end 

(21) 


which,  together  with  (9).  relates  the  time  dependence  of  the  dynamic  stress 
intensity  factors  using  deformation  theory  co  ch*  dm*  dependence  of  Che 
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corresponding  intensity  factors  using'  linear  elasticity. 
To  understand  the  implications  of  Eqn.(21),  let  f(t) 

'P(a+l)  - - 

_N6 


f(D) 
for  f 


r(crt-l) 


P 

t.f5  «e  ft(t)  - C c*.  yt>  ■ 
form  of  f(,  ij  into  Eqn.(2l)  gives 


where 

A 


is  the  g< 

,N 


n+i  r(6N+n  m+n 

5N+1  5+1 

p p 


F5H 


c*  (a  > 0).  i.e. 

function.  The  corresponding  forms 
Substituting  the  Laplace  trans- 

(22) 


from  which  we  find  5 • 2cr/(N+l).  Since  0 < S s 1 , it  follows  that  5 > a 
which  implies  that  ft(c)  ~ increases  slower  than  f(t),  but  f (t)  ~ t* 
increases  faster  than  for  the  linearly  elastic  case. 

INFLUENCE  OF  CONSTITUTIVE  BEHAVIOR  ON  NEAR-TIP  FIELDS 


The  functional  dependence  on  r and  9 of  the  fields  of  stress  and  de- 
formation near  a rapidly  propagating  crack  tip  can  generally  be  established 
by  asymptotic  considerations.  For  a wide  claas  of  constitutive  behaviors  an 
asymptotic  analysis  can  be  based  on  "separation  of  variables"  near-tip 
solutions,  in  which  the  dependence  on  r is  a-priori  assumed  as  a power  of  r. 
The  analysis  proceeds  by  collecting  the  most  singular  terms  in  the  governing 
equations,  and  the  boundary  conditions,  and  subsequently  solving  the  resulting 
linear  or  nonlinear  eigenvalue  problem.  When  this  method  is  applicable  It  is 
found  chat  the  stresses  and  strains  are  singular. 

In  an  asymptotic  analysis  only  the  field  in  the  highly  strained  material 
in  the  near-tip  region  is  considered.  Hence,  it  is  the  stress-strain 
behavior  at  very  large  strains  which  enters  an  asymptotic  analysis.  The 
speed  of  crack  propagation  is  limited  by  che  characteristic  speed  of  the 
material  insedlately  ahead  of  the  crack  tip.  This  speed  is  related  to  che 
slope  of  the  stress-scrain  curve,  Thus , only  constitutive  models  with  a 
finite  slope  of  the  stress-strain  curve  at  large  scrains  are  suitable  for 
the  type  of  asymptotic  analysis  described  above. 

In  this  section  Che  Influence  of  che  constitutive  model  on  che  near- dp 
field  is  examined  for  rapidly  propagating  cracks.  Only  the  case  of  Mode  III 
crack  propagation  (antiplane  strain)  has  so  far  been  investigated  in  some 
decail. 

The  asymptotic  analysis  is  carried  out  in  a moving  coordinate  system, 
which  is  fixed  with  respect  to  che  crack  tip.  If  che  speed  of  crack  propa- 
gation is  v » dZ(c)dt,  where  Z(C)  is  crack  length  and  v is  an  arbic-ary 
function  of  time  subject  to  the  conditions  chat  v and  dv/dt  be  continuous, 
the  material  time  derivatives  are  transformed  to  Che  moving  coordinate 
systems  by  the  relations 


( ) 


if  - v(t) 


and 


••  a >2 

< ) ■ 77  - *<e)  f?  * 2 v(c)  Sefe 


♦ v2(t)  -5 

1 *X2 


(23) 


In  the  usual  separation  of  variable  type  asymptotic  analysis  considered  here, 
the  terms  coming  from  u./5xj  will  be  more  singular  than  any  ocher  terms 
in  ii  , so  that  the  results  will  be  the  same  for  the  "steady  state"  and 
transient  cases  as  far  as  che  r and  9 dependence  is  concerned.  By  "steady 
state",  it  is  meanc  chat,  as  seen  by  a crack-tip  observer,  the  stress  and 
strain  fields  are  constant. 

Motion  in  antiplane  strain  is  defined  by  a displacement  in  the  x«- 
direction  only,  see  Fig.  1.  The  notation  is  simplified  by  writing, 
w ■ u (x, ,x„)  , u » u ■ 0 , 

Xj  12  Xj  Xj 


6 


1 - 1,  2 


T • C 

1 31  ’ 

Vt  - 2«31  • . 1 • 1.  2 (24) 

The  equation*  of  motion  and  coaipacabillty  reduce  to 

Ti,i  " 0 " > Yl,2  ' y2,1  1 (25) 

where  w Is  found  from  Eqn.(23). 

The  system  of  governing  equations  must  be  completed  by  the  constitutive 
relations.  In  the  following  we  examine  the  near- tip  fields  for  the  various 
different  constitutive  behaviors  shown  In  Fig.  2. 


Fig.  2 Models  of  material  behavior 


Linear  Elasticity  (Fig.  2a) 

The  constitutive  law  Is  r » l*  y and  the  near- dp  field  Is  of  che 
fora  given  in  section  2.  Thetis  the  stresses  and  strains  have  the  familiar 
square  root  singularity  with  respect  to  r,  of  LEFM.  The  upper  limit  for 
crack  propagation  under  remotely  applied  loads  is  the  elastic  shear  wave 
speed  <>/o)’. 


Bilinear  elasticity  (Fig.  2b) 

Defining  an  affective  shear  stress,  t,  and  shear  strain,  v,  as 
T » (r^  T^)v  i Y ■ (y^Y^)*  , the  constitutive  law  is  given  by 


kk  « Tt  + ( n/nc  - D(1  - t0/t)  Tt  , T > ro  (26) 

where  tk  is  che  linear  elastic  shear  modulus  and  a * dc/dv  for  t > - 

C 0 

The  singularities  for  stress  and  strain  are  similar  to  the  ones  for 
classical  linear  elasticity,  but  with  u replacing  |k  as  che  relevant  elastic 
constant.  This  reduction  in  the  slope  of  che  stress-strain  curve  has 
important  consequences  for  che  significance  of  dynamic  effects  in  rapid  crack 
propagation.  Although  the  crack  speed  v may  be  small  compared  to  che  linear 
elaaclc  shear-wave  speed,  (,k/c)  , it  may  be  a significant  fraction  of  the 
characteristic  wave  speed  at  high  strains,  (tkt/c)%and  this  is  the  important 
comparison  in  deciding  the  importance  of  dynamic  effects.  If  the  fracture 
process  is  essentially  brittle,  (i*  /i)'  is  an  approximate  upper  limit  on  che 
crack  propagation  velocity. 


Power-law  Strain  Hardening  (Fig.  2c) 


i 


The  Mod a III  version  for  deformation  theory  of  plasticity  Is 


T1  ' * Yi  ’ 


T « T 


Ti  “ T0  ( Y/Yo  r Yi  ’ T > To 


(27) 


where  r , are  the  yield  stress  and  strain  respectively  and  H > 0. 

Die  loading  path  In  stress  space  of  a material  point  close  to  the  tip  of 
a rapidly  propagating  crack  la  not  proportional  since  a sone  of  unloading 
exists.  Deformation  theory  of  plasticity  will  therefore  be  a poor  model  in 
this  region.  Even  If  Eqn.(27)  is  assumed  to  be  velid  for  a non-linear  elastic 
material,  it  is  not  possible  to  obtain  an  asymptotic  expansion  in  potters  of 
r.  For  0 < N < l , ii  (t)  • dr/dy  vanishes  is  y * ■ , and  the  characteristic 

wave  speed  becomes  zero.  Therefore,  any  speed  of  crack  propagation  Is  super- 
sonic and  the  above  asymptotic  expansion  is  not  valid.  For  S > l , dr'dy 
becomes  unbounded  as  y “ » and  so  does  the  characteristic  weve  speed.  This 
means  that  for  a bounded  crack  tip  velocity,  dynamic  effects  disappear 
altogether  for  N > 1. 


J2  Flow  Theory  (Fig.  2d) 

Finally,  a crack  propagating  rapidly  in  a strain-hardening,  elastic- 
plastic  material  characterized  by  J2  flow  theory  and  a bilinear  stress- 
strain  curve  is  considered.  For  loading  in  the  plastic  regicrs,  Che 
incremental  s cress- strain  relations  are 

*c  Y1  " * fl  + (l  - Ti  ’ T > To  " + > 0 ’ (28) 

while  for  elastic  unloading  and  elastic  loading  we  have 

» y « 0 + , TiO  , or  TST  , (29) 

c 1 1 o 


where  a * * '*  • 

If  the  sceady- state  situation  is  considered,  the  ebove  equations  toget- 
her with  the  equations  of  mocion  and  compataoility  arc  not  alvays  elliptic. 
When  the  governing  equations  are  elllpdc,  an  asymptotic  expansion  in  powers 
of  r will  be  valid.  If  the  material  time  derivatives  are  replaced  by 
• v 5/5*1  and  the  comparability  equation  in  Eqn.(25)  is  used  to  eliminate  vj 
from  the  constitutive  equation  (28),  the  governing  equations  in  Che  plastic 
loading  region  can  be  written  as  Che  following  system  of  first  order 
equations , 


l 

- U + (1/a  - l)]/(pT2) 

- T1T2(l/®  * D/(*t2) 


TjTjd/q  - l)/(*0 
[l  + r2(l/a  - 1)]/(*t2) 


where 


+ 


0 

0 

1 


0 

0 

0 


(30) 
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Using  the  fora  A + A a),«  ■ 0 to  define  the  matrices  A , A , the  charac- 
teristic directions  l[l_3]  (where  characteristics  are  defined  as  lines  across 
uhlch  discontinuities  in  j>  can  occur)  are  given  by  A • dxj/dx^.  these 
directions  are  found  by  solving  the  characteristic  equation 

|A2  - A A1  | • 0 . (31) 

the  system  of  equations  is  elliptic  if  the  roots  A are  complex. 

Equation  (31)  can  be  written  in  the  fora  Af(A)  • 0,  where  f(A)  is 
quadratic  in  A,  with  roots  A~  which  will  be  complex  if  the  discriminant  of 
f(A)  is  negative,  this  condition  can  be  reduced  to 

[v2/(ii/o)]Cl  + t2  (1/a  - l)/t2]  - 1 < 0 (32) 

If  “ T , a 'ower  bound  of 

v < (*t/p)*  (33) 

is  found  for  che  above  system  of  equations  to  be  elliptic  ■ 14 J.  which  is 
the  limit  given  earlier  on  the  basis  of  intuitive  reasoning. 

For  crack  dp  velocities  less  chan  this  limit,  the  separation  of  varia- 
bles approach  to  an  asymptotic  expansion  will  be  valid.  However,  for  higher 
velocities,  the  type  of  the  governing  equations  depends  on  the  magnitude  of 

T1  • 1. 

the  case  of  v < (u  /o)  will  be  considered  below,  and  in  che  next  section 
che  elastic  perfectly-plastic  case  will  be  considered.  Following  Achenbach 
and  Kanninen  [6],  an  asymptotically  valid  solution  for  w of  the  general  fora 

w - C v W(6)  rs  (34) 

is  sought.  Here  C is  an  amplitude  factor,  left  undetermined  by  the  analysis, 
and  K(9)  and  s have  to  be  found,  the  problem  is  set  up  as  a generalized 
eigenvalue  problem  with  s as  the  eigenvalue. 

Using  Eqn.(24),  the  strain  rates  corresponding  to  Eqn.(33)  are 

vt  • c V [w(9)  r*]  . i - 1,2  (35) 

Since  only  the  mosc  singular  terms  in  r will  be  retained  in  che  asymptotic 
analysis,  in  Eqn.(23)  quantities  of  the  fora  3f/>t,  where  f is  any  field 
variable,  may  be  neglected  when  compared  with  -v(t)  if/Sx,,  so  that  ( ')  m 
- v c/9x^.  In  the  asymptotic  analysis  v is  considered  to  be  constant, 
the  stresses  corresponding  to  Eqn.(24)  are 


where  superscript  ' indicates  3/33  and  32  ■ y2/(i*/p). 

Substitution  of  Eqns.(33)  and  (37-39)  into  Eqn.(28)  yields 


a(s  H cos9  - W slnd)  ■ » Tj  + (1  • j)  I Ij  I 


-1 


s(s  W sine  + w cos0)  » a t2  + (1  - a)  T T2  T 


(41) 

(42) 


The  corresponding  equations  for  elastic  unloading  can  be  obtained  from  [15]. 
The  pertinent  s~’ution  is 


• O *>  m / 2 

W • (1-3  sin‘3)  ' cos  [s(ui-n)]  , 


(43) 


where  tan  tu  » (1-3  )’  tan  9 


The  problem  is  and- symmetric  with  respect  to  9 • 0 , so  only  the  domain 
0 s 8 s n Is  considered.  The  boundary  conditions  are  given  by 


W 

W' 


on  9 » 0 (anti-symmetry)  , 
on  8 ■ it  (zero  tractions) 


(44) 


As  in  the  quasi-static  solution  \ 16 'i  , the  possibility  of  a plastic  reloading 
zone  along  the  crack  flanks  is  not  considered.  Although  this  may  be  a 
reasonable  assusiptioa  for  Mode  III,  results  for  the  quasi-static  elastic 
parfectly-plastic  Mode  t case  [ 17 1 indicate  that  the  above  assumption  is  not 
correct  in  Mode  I. 


Fig.  3 Boundary  between  loading  and  unloading 
regions  as  r - 0. 


The  boundary  between  loading  and  unloading  zones  surrounding  the  crack 
tip  is  assumed  to  be  a radial  line,  emanating  from  the  crack  tip  at  an  angle 
5 » 9 , (Fig.  3).  The  field  in  the  loading  zone,  0 s 3 s 9 , is  governed  by 
Eqns.p(40)  and  (41-42)  and  the  field  in  the  elastic  unloading  zone 
c s 3 s rr  is  given  by  Eqn.(43). 


To  completely  specify  the  problem,  the  continuity  conditions  at  the 
Interface,  3 » 3p,  must  be  given.  From  the  unloading  condition,  t » 0 , 
which  implies  T » 0 , we  find 


- s T cos6  + T'  sin0  » 0 at  3 


(45) 


In  addition,  it  is  assumed  that  the  particle  velocity  and  stresses  are 
continuous  at  3 » 3 . This  condition  can  be  written  as 
P 


[W] 


at  9 • 3 


(46) 


where 


The  first  condition  of  Eqn.(46)  can  be  automatically  satisfied  by  writing  the 
solution  in  the  elastic  unloading  zone  as 


H,  * W(9‘>  M,(9)  / W,(8  ) (47) 

where  17(3^)  is  the  solution  in  the  loading  region  at  3 • 3^.  From  continuity 

of  W at  3 ■ 3 , it  follows  that 

P 

H'(8')  - W(e')  W'(8)  / W (8  ) , (48) 

P P • P « P 

where  W (8  ) and  M'(8  ) are  found  from  Eqn.(43). 

• P e p 

The  problem  has  now  been  reduced  to  determining  the  solution  in  the 
plaatic  loading  region,  l.e.  a solution  satisfying  the  field  equations  in  the 
region  0 s 3 s 3 , and  the  boundary  conditions  at  3 • 0 , 6«  . The 

quantities  whicti  have  to  be  determined  are  W(8),  Ti(6)  , T?(8),  3.  and  s. 

This  is  a nonlinear  eigenvalue  problem  which  must  be  solved  numerically. 

The  results  for  s and  3 are  shown  in  Tables  1 and  2,  and  further  details  of 
the  analysis  are  given  In  [6]. 

The  limit  on  v derived  earlier,  is  an  approximate  agreement  with  the 
region  where  the  numerical  solution  of  the  equations  given  above  could  be 
found. 


TABLE  1.  CALCULATED  VALUES  OF  3-  FOR  DYNAMIC 

PLASTIC  ANTI- PLANE  SHEAR  CRACK  PROPAGATION 


3 

3-0 

3-0.1 

3-0.25 

3-0.5 

3-0.75 

1.0 

1.571 

1.576 

1.602 

1.690 

1.786 

0.7 

1.522 

1.528 

1.554 

1.643 

1.731 

0.5 

1.473 

1.478 

1.505 

1.595 

0.3 

1.393 

1.398 

1.427 

1.519 

0.2 

1.328 

1.334 

1.363 

0.1 

1.217 

1.225 

1.259 

I 


TABLE  2. 

CALCULATED 
ANTI -PLANE 

VALUES  OF  s 
SHEAR  CRACK 

FOR  DYNAMIC 
PROPAGATION 

PLASTIC 

t 

or 

3-0 

3-0.1 

5-0.25 

3-0.5 

3-0.75 

1.0 

-0.500 

-0.500 

-0.500 

-0.500 

-0.500 

0.7 

-0.444 

•0.444 

-0.442 

-0.434 

-0.396 

0.5 

-0.395 

-0.394 

-0.391 

-0.375 

0.3 

-0.325 

•0.324 

-0.319 

•0.2S8 

0.2 

-0.277 

-0.276 

-0.269 

0.1 

•0.208 

-0.206 

-0.194 

MODE  III  STEADY-STATE  DYNAMIC  CRACK  PROPAGATION 
IN  AN  ELASTIC  PERFECTLY -PLASTIC  MATERIAL 


la  the  previous  section,  the  case  where  the  system  of  governing  equations 
is  elliptic,  and  a sepsrstion-of- variable  type  asyoptotic  analysis  is  appli- 
cable, was  considered.  In  this  section  the  special  case  of  an  elastic 
perfectly-plastic  material,  for  which  the  governing  equations  in  the  plastic 
loading  region  are  hyperbolic,  is  investigated.  The  analysis  is  again 
asymptotic  in  the  sense  thet  the  solution  is  not  velld  far  from  the  crack  tip. 
However  it  is  different  from  the  approach  where  only  terms  of  the  highest 
singularity  in  r are  retained. 

Ihe  problem  is  a generalisation  to  dynamics  of  the  problem  addressed  by 
Chitaley  and  McClintock  [£]  for  quasi-static  crack  growth.  The  Mode  III  case 
is  used  as  a means  to  understand  the  importance  of  inertial  terms  on  the 
crack  tip  fields,  and  to  gain  insight  into  the  more  difficult  Mode  I case. 
Although  for  an  analysis  of  a real  macerlal,  the  exclusion  of  strain  rate 
effects  in  the  material  model  cannot  be  Justified  when  other  dynamic 
plasticity  effects  are  Included,  this  simplification  is  made  to  make  the 
problem  tractable  to  an  asymptotic  analysis.  The  results  which  are  obtained 
on  the  basis  of  several  assumptions  do  show  that  when  compared  to  the  quasi- 
static  case  [9],  there  is  a major  difference.  As  opposed  to  the  cases  dis- 
cussed in  the  previous  section  it  is  not  clear  that  the  equations  governing 
the  asymptotic  field  for  the  transient  case  are  the  same  as  those  for  the 
steady-state  case  when  the  material  behavior  is  elastic  perfectly-plastic; 
chat  is,  it  cannot  be  assumed  that  ( " ) ~ — v 3/3x^  is  a good  approximation 
for  the  transient  case.  However,  only  the  steady-state  case  will  be  con- 
sidered here. 


Constitutive  Law  (Fig.  2d) 


The  yield  condition  is 
2 

T * 

o 


(49) 


Adopting  a Prandtl-Reuss  incremental  flow  law  as  in  Eqn.  (26),  the  strain 
rates  are  given  by 


V*  + X Ti 


(50) 


where  X.  > 0 if  r * r , r > 0 , otherwise  X * 0 . Equation  (50)  Is  a special 
case  of  Eqn.(28)  but°as  the  analysis  Is  fundamentally  different,  the  governing 
aquations  ara  radar ivad. 

Governing  Equations 


Using  the  steady  state  form  of  Eqn.(23),  i.a.  all  time  derivatives,  ara 
assumed  to  vanish  Identically,  the  governing  equations  can  be  written  In  the 
following  form. 


i - p v2  32w/5x2  , 

(equilibrium) 

(51) 

Y1.2  " y2,1 

(compatabllity) 

(52) 

Yi,l  “ Tl,l/l*  + X’l  Ti 

(constitutive  law) 

(53) 

where  X,.  < 0 If  r « t , t,  t,  , 
l o 1 1,1 

< 0 , 

otherwise  X>2  » 0 . 

In  the  plastically  deforming  regions,  where  X,  <0,  the  yield  condition 
can  be  Identically  satisfied  If  we  define  ® such  that 


T.  - - T sin® 
l o 


T COS# 

o 


(54) 


Substituting  into  Eqns. (51-53)  and  eliminating  y,  , the  governing  equations 
can  be  written  as 


T 

COS® 

X* 


T • T„ 

— cos®  — r 

**  *2 


*’l  + 


— sin® 

W 


sin® 


£•2*2;  a>  “ , 


(55) 


1 2 

where  3 “ v/^Wp  . Equation  (35)  has  the  general  form  A + A m.j  ■ 0 • 
Using  the  method  of  characteristics  .13],  Eqn.(55)  Is  transformed  Into 
a set  of  ordinary  differential  equations  along  the  characteristics  In  the 
x,  - x»  plane.  To  do  this  Is  written  as  • daj/dxj-ldxj/dx, ) m,, 
where  dm/dx|  indicates  differentiation  along  the  characteristic  whose  tangent 
is  given  by  dxj/'dxj.  Substituting  this  relation  into  Eqn.  (55)  results  In 

A1  d£/dx1  + (A2  - A1  dXj/dXj)  m>2  - 0 . (56) 

2 1 

Defining  A as  A.  • dx2/dx^  and  solving  for  A from  |A  - A A I • 0 , which  In 
this  case  Is  given  by 

T2  T2  T2 

A2  (22-cos2®)  + 2 A -|  sin2®  cos2®  - sin2®  • 0 , (57) 

* U It 


the  characteristic  directions  are  found  to  be 


lis2_ 

cos®  s 2 ’ 


(58) 


Here  the  - signs  are  used  to  denote  the  = ve  characteristic  curves.  The  left 
eigenvectors corresponding  to  A,whlch  ara  defined  by  ] A*  - A A1]  • 0 , 

are 
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If  1 Is  contracted  vlch  Eqn.(55),  th«  result  Is  sn  ordinary  differencial 
equation,  called  che  characteristic  relation,  along  each  characteristic, 
wit, 

l. A1  dui/d*1  + X.(A2  - .'A1)  a,2  • i.A1  duj/dx1  - 0 . (60 

For  che  problem  at  hand  the  characteristic  relations  are 

6 dy.  s (T  /»)  d®  - 0 (61 

L o 

Equation  (61)  can  be  integrated  along  the  charecterlstlcs  to  give  the 
Rlemann  invariants, 

J*  ■ = (Ww)  « (62 


where  are  constants  along  the  corresponding  characteristics. 

From  Eqns.(61)  and  (62)  a number  of  cypes  of  characteristic  fields  in  a 
region  in  che  Xj  - X2  plane  can  be  constructed.  They  are 

a)  simple  field  , J+  - constant,  ■ constant, 
both  sets  of  characteristics  straight; 

b)  uniform  field,  J+  • constant,  J_  4 constant, 

-ve  set  of  charecterlstlcs  straight; 

c)  uniform  field,  4 constant,  • constant, 

+ve  set  of  characteristics  straight; 

and 

d)  non-uniform  field,  J+  4 constant,  J_  4 constant, 
neither  set  of  characteristics  straight. 

At  this  polnc,  che  difference  between  che  governing  equations  in  the 
plastic  zones  for  che  dynamic  and  quasi- static  cases  can  be  pointed  out.  If 
in  Eqns.(61)  and  (62)  che  limit  3 - 0 is  taken,  the  two  characteristics  and 
the  corresponding  Rlemann  invariants  dagenerate  to  one  charscterlstlc  and  or.e 
Rlemann  invariant,  and  the  appropriate  quaai-statlc  limit  is  obtained.  The 
resulting  set  of  equations  is^lles  chat  ® is  constant  along  a characteristic 
and  therefore  che  characteristics  are  straight.  Unfortunately  in  the 
dynamic  case  under  consideration,  che  characteristics  are  in  general  not 
straight  since  their  direction  at  any  point  depends  on  the  solution,  0 , at 
chat  point. 


Boundary  Conditions 

The  Mode  III  problem  is  antl-synsetrlc  ebout  xj  " 0 , which  implies 
w • 0 along  x,  ■ 0 , x^  > 0.  Therefore  Yj  • 0 on  the  X}  - axis  ahead  of  the 
crack  and  by  Eqn.(53),  this  implies  ■ 0 or  3 ■ 0 on  Xj  • 0 , x^  > 0 
On  the  crack  flanks  (x<  ■ 0 , x^  < 0)  we  have  Tj  » 0 . The  boundary 
conditions  near  the  crack  tip  are  then 

’•2  * < 


2 


In  Eqn.(64),  eh*  choice  of  being  tq  or  -t^  depend*  on  the  solution 


given  below. 

The  boundary  conditions  far  froa  the  crack  tip  are  not  specified  as  they 
will  not  enter  the  analysis.  However,  It  is  envisioned  chat  the  snail  scale 
yielding  assumption  will  hold  and  the  plastic  zones  will  be  embedded  In  the 
linear  elastic  stress  field  for  a moving  crack. 

Since  the  problem  Is  anti-symmetric.  only  the  upper-half  plane, 

X2  > 0 , will  be  considered. 


Hear-Tip  Fields 


From  Eqn.(62),  a difference  can  be  seen  between  the  quasi-static  and 
dynamic  cases.  From  19,1^,17],  the  quasi-static  scraln  y^  is  found  to  be 
logarithmically  singular  as  the  crack  tip  is  approached.  However.  In  the 
dynamic  case.  If  yi  should  become  unbounded  as  r - 0 but  4 is  restricted  to 
finite  values.  It  follows  from  Eqn.(62)  that  oust  also  tend  to  infinity. 
The  only  way  this  could  happen  with  bounded  4 Is  for  4 to  be  cyclic  in  its 
permissible  rang*.  The  latter  would  imply  that  tha  stresses  oust  also  be 
cyclic,  which  is  unacceptable  for  a solution  in  a plastic  zone  Thus,  is 
assumed  to  be  bounded  la  the  dynamic  css*. 

A solution  is  now  constructed, which  satisfies  the  boundary  conditions  on 
*2  “ O.by  using  Che  simple  and  uniform  fields  a)  and  b)  given  above. 

Although  it  is  possible  that  the  near-tip  plastic  loading  field  may  be  non- 
uniform,  this  case  is  not  considered.  Case  c)  is  ruled  out  since 
V-i  > 0 and  cases  a)  or  b)  taken  separately  cannot  match  both  boundary 
conditions.  Tha  case  of  an  elastic  unloading  zone  is  also  considered,  but 
will  be  shown  not  to  exist  to  the  order  of  approximation  made. 

After  a process  of  elimination  the  following  solution  is  found.  The 
resulting  characteristics  are  shown  in  Fig.  4. 


Fig. 


4 Characteristic  curves  near  the  crack  tip 

for  Mod*  III  steady-state  dynamic  crack  propagation 


Let  a simple  field  border  on  Xg  * 0 , „ - _ 
characteristics  are  straight,  being  given  by 


xt  < 0. 


Then,  in  this  region  the 


= Xj/3  + Ax 


(65) 


where  A±  are  constants  and  0 , V|  are  constant.  Applying  the  traction  free 
crack  face  boundary  condition  given  in  Eqn.(64),  results  in  4 ■ ">'2  in  the 
simple  field.  The  choice  of  4 > 0 will  be  obvious  once  the  field  in  front 
of  tha  crack  tip  is  found . 

y^  • 0 satisfies  the 
axis  is  a charac- 


Tn  front  of  the  crack  tip  , on  the  x^  - axis  , 4 


characteristic  relations  in  Eon.  (61).  Therefore  the  xj  - axis  is  a 
terlsclc.  Also,  from  the  stress  distribution  in  Che  linear  elastic  case, 
a plastic  loading  zone  is  expected  to  form  ahead  of  the  crack  tip.  Using 
this  as  motivation,  case  b),  with  the  -ve  characteristics  straight  and 
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forming  a can  trad  fan  at  eha  crack  dp  ia  a a ad  to  match  cha  boundary  con- 
dition ahaad  of  cha  crack  tip.  Caaa  b)  aaauaaa  J+  ia  conacane,  which  ra- 
quiraa  Y],  , 3 be  conacane  along  a -va  characteristic . Substituting  this  into 
Eqn. (58) , cha  -va  characeariacics  ara  found  co  be  defined  by 


x2  ' X1  coed  - 3 


x2  * xl 


where  cha  conscant  of  integration  la  taro  to  sadafy  eha  eencrad  fan  con- 
dition, and  3 is  cha  angle  bacvaan  the  characteristic  and  the  positive 
- axis. 

Differentiating  Eqn. (66)  with  respect  co  x^  and  solving  for  3 in  earma 
of  3 from  eha  raaulting  aquation,  wa  find 

coa  3 • 3 sin23  + cos3  (1  - 32sin29)*  . (67) 

Substituting  this  result  into  Eqns.(54)  gives 

rx  - -tq  sing  [(1  - 32  sin29)'  - 3 cos9]  , (68) 


[3  sin29  + cos8  (l  - 32  sin28)^]  . 


In  the  simple  field,  Tj  ■ 0.  The  angle  at  which  the  uniform  and  simple  fields 
match  is  found  from  setting  Eqn. (69)  to  taro.  The  result  is 

can8*  - - 1/3  , (70) 

which  defines  a characteristic  of  both  the  simple  and  uniform  fields  emanating 
from  the  crack  dp. 

Since  the  x, - axis  is  a characteristic  on  which  v.  • ; j 0 , J+  lj  zero. 
Therefore,  from  Eqn.  (62),  Vj“  - to®/(I*6)  and  using  Eqn.  (67),  we  can  write 


vl  * - jjf  cos  2 [3  sln29  + cos8  (1  - 32  sin23)^]  , 

* 

for  0 £ $ £ 8 , 

★ 

wh«r«  d Is  defined  by  Eqn.  (70).  From  Eqn.  (53) , a.,^  cnn  b«  obtained  as 

x . COS^)  Hi 
’1  t r d8  ’ 
o 


(71) 


(72) 


and  by  virtue  of  Eqn. (71)  it  can^then  be  easily  shown  that  < 0 as 
required.  By  evaluating  v,  at  3*.  its  value  in  the  simple  field  can  be 
obcalned  as  • (-T  roy(2ii=).  If  Eqn.  (70)  is  substituted  in  Eqn.  (68),  it  is 

found  chat  r.  • -t  , 0 " t/2,  as  previously  assumed. 

1 o 


The  possibility  of  an  elastic  unloading  sons  extending  to  the  crack  dp 
between  the  two  plastic  fields  must  now  be  checked.  If  the  material  unloads 
such  chat  V,i  > 0 , chere  will  be  a wake  region  with  residual  plastic  strains 
behind  the  plastic  loading  zone.  Obviously  these  residual  plastic  strains 
will  only  be  functions  of  The  governing  equations  in  the  wake  region  are 


where  Yj  la  cha  non- recoverab La  pare  of  tha  y ^ aeralo.  F 
evaluated  ae  6 • 8*  , le  can  ba  aaan  that  y2  and  '■  ae  6* 


From  Eqn. (71) . 
ara  not 

function*  of  x2  and  eharafora  In  cha  elastic  unloading  raglon  2 la  aaro. 

No  aolutlon  can  chan  ba  found  of  Eqns. (73-74)  which  natch**  cha  continuity 
condition*  on  both  aldaa  of  eh*  alaaelc  unloading  raglon,  and  eharafora  cha 
assumption  of  an  unloading  aona  la  lncorract. 

From  Eqn.(S2),  Y2,i  can  b|  found  and  afear  lneagraelng  with  raapace  to 
x,,  y2  In  tha  donaln  0’s  i s J la  obeainad  aa 

Yj  ■ In  [l  - 3 sln29  - cos8  (l  - 32sin29)^] 

+ •—  In  [l  + 3 aln29  + cos8  (i-32sin28)'] } + s(x2>  . (73) 

whar*  -P(x2)  will  b*  apaclflad  laear.  By  continuity  of  aeralna,  for 

9*  S 8 s tt 

v2  - <p(x2)  (76) 

Requiring  Yj  to  ba  bounded  at  x1  • r and  for  small  Xj.  where  x^  • r la  the 
lntaraacelon  of  cha  alaaelc-plaatlc  boundary  with  the  X|-  axis,  «(x2)  can  b* 
approximately  given  by 


,(X2)  " M la  W * Xn(l*3)  * in2  + l} 


(77) 


by 


From  eh*  y2  a train  In  the  a Imp la  field,  the  crack  opening  angle  la  given 


9 « 2 can 

op 


-1 


,Tt  T 
\ 2w3/ 


(78) 


and  the  crack  opening  dlaplacamane  la  aaro.  If  r la  cha  radial  dlatanc*  co 
a point  on  a -Hr*  charactarlatlc,  the  -Hr*  characteristic*  are  given  by 


dr  . , (3  co«9  ,-t.  .a-,32si.n23)\* 
r \ 23*ln6  / 


,79) 


Tha  above  result*  ara  aaaaneially  the  same  aa  thoa*  found  by  Slapyan 
.18.,  who  uaad  an  approach  whar*  only  cha  most  singular  carma  (0(l/r))  w*ra 
retained  In  tha  governing  aquations.  This  approach  can  b*  shown  to  glv*  the 
same  result  as  whan  cha  field  variable  are  considered  to  b*  functions  of 
(x^/xj)  only.  Slnca  tha  results  given  above  agree  with  ll§],  they  muse  also 
b*  solutions  for  cha  case  whar*  the  only  Indapandane  variable  la  (x^/xj) ■ 

Ic  la  Intuitively  reasonable  chat  eh*  steady-scat*  dynamic  results 
should  reduce  Co  Cha  quasi-static  results  ae  small  values  of  3.  However  the 
cransicion  from  dynamic  to  quasl-statlc  does  noc  seem  to  be  uniform  slnca  If 
3 - 0,  Yi  and  Y2  Increase  beyond  bounds. 

In  the  quasl-statlc  case  .9,17],  y2  ~ 2n(r)slnd  as  r - 0 and  char*  Is 
a raglon  of  elastic  unloading.  As  "a  material  point  passes  tha  crack  tip,  eh* 
Increments  of  v2  changa  sign  and  this  can  only  b*  achieved  by  elastic  un- 
loading in  the  quasl-statlc  casa.  The  same  affect  Is  not  found  In  tha 
dynamic  casa,  slnca  the  plastic  loading  region  ax tends  completely  around  eh* 
crack  tip.  An  estimation  of  the  relative  son*  sis*  near  the  crack  tip  over 
which  the  dynamic  results  ara  valid  can  ba  found  from  Eqns. (35) . Sire* 
inertial  affects  ara  proportional  to  3*  , the  dynamic  results  will  ba 
Important  In  a raglon  of  0(3^)  , Indicated  In  Fig.  4 by  Cha  minimum  distance 
between  the  speculated  clastic-plastic  boundary  and  tha  crack  tip.  Outside 


Tabla  3: 


Comparison  of  results  from  Ref . [l2,9, 19]  and  the  praaanc  analysis 


I 


I 


seaclonary 

crack 

(static) 


propagating  crack 
(staady-stata) 


quasi-static 


dynamic 


acraln  y2 

- - 0,  Xj  > 0 


Sjft) 

Raf. [l2] 


Yo[1+ln(~^) + i /q2  H)j 

Ref . [l2] 


*■&» 


u^)! 

Raf. [12] 


-i  (W 

o 

Raf. 19] 


l/^II)2  l 

TT\  T 


1 (1-22)* 


crack* opening 
angle  8 

op 


2 e“'l(-5 

Raf. [12] 


2 tan 


Raf.  [19] 


crack* opening 
diaplacamant 


2 <*?!!>* 
"WT- 


Raf. [12] 


Raf.  [19] 


S 

■ static  or  quasi-static  stress  intensity  factor 


"ill 


■ dynamic  stress  intensity  factor 


A comparison  is  made  in  Tabla  3 between  the  static,  quasi*statlc  and 
iynamlc  results  for  a semi-infinite  crack  and  small  scale  yielding.  The 
approximate  plastic  zone  size,  r_,  for  the  dynamic  case  is  estlaiated  by 
using  the  Irvin  approach  of  an  effective  crack  length  [20]  and  does  appear 
to  be  consistent  with  the  quasi-static  result,  although  this  may  be  purely 
fortuitous.  In  the  limit  as  3 - 1,  k!^  - 0 so  that  rp  - 0 which  would  be 

the  correct  limit  for  the  supersonic  case.  The  crack-opening  angles  and 
displacements  are  consistent  with  e00  - n as  2 - 0.  The  Y2  strains  along  the 
x-axis  are  not  consistent,  but  this  is  not  unexpected,  as  the  region  over 
which  Che  dynamic  solution  dominates,  tends  to  zero  in  the  limit  2-0.  It 
is  expected  that  the  strain  outside  the  region  in  which  the  present  analysis 
is  valid  will  tend  to  the  quasi-static  result. 

The  above  analysis  has  taken  no  account  of  loading  or  elastlc-plaatic 
boundary  effects  aod  is  based  on  a number  of  simplifying  assumptions.  By  a 
process  of  elimination,  the  above  solution  has  been  found  to  be  Che  only  one 
satisfying  these  assusipcions . Further  detailed  analysis  requires  the  help 
of  extensive  numerical  work  to  give  guidance  in  the  analytical  work. 
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